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Abstract 


This final report summarizes the research carried out on the development of numerical 
approaches for simulating the radiation emitted in high-altitude plumes of spacecraft and 
missiles. In Appendix I, we describe the development of ultra-violet radiation model for the 
hydroxyl radical that is implememted in a particle-based Monte Carlo approach for simulating 
the plumes. In Appendix II, we describe a plume radiation calculation for the plume of the 
ammonia arcjet flown on the USAF ESEX flight experiment. The research performed in this 
grant has led to significant advances in the capabilities for modeling these flows. 
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Overview 


It was the aim of this research to improve prediction of ultra-violet and visible radiation emitted 
by high altitude missile plumes and by the interaction of spacecraft control thruster plumes with 
the ambient atmosphere. To improve prediction of radiative emission it is first necessary to 
establish accurate modeling of the fluid dynamics and chemistry for these flows. This research 
employs the direct simulation Monte Carlo method (DSMC) to compute the flows that are in a 
strong state of thermochemical nonequilibrium. For application to these flows, development of 
the DSMC technique is required in terms of electronic excitation models and in terms of schemes 
for including trace species. These issues were addressed in this research program. 

In the course of the grant, significant progress was made in two areas. The first concerned the 
development of DSMC for direct prediction of ultra-violet radiation. Technical details are 
provided in Appendix A which contains a reprint of a paper published in the Journal of 
Spacecraft and. Rockets. Briefly, the main accomplishments were: (1) direct simulation of the 
chemistry of formation of the hydroxyl excited state OH(A); (2) detailed prediction of ultra¬ 
violet radiation from hydroxyl for hypersonic flows; (3) direct verification of the kinetics models 
against experimental data. 


The second area of work was conducted in support of the USAF flight experiment called ESEX 
(Electric propulsion Space Experiment). On ESEX, a high power ammonia arcjet was fired in 
space. Flight measurements on the spacecraft included contamination and thermal radiation. 
Remote observation of the radiation emitted by the arcjet plume were made from a USAF facility 
in Maui. Detailed DSMC computations were performed using 10 million particles on 32 
processors of an IBM SP-2 parallel supercomputer in order to accurately resolve the spacecraft 
contamination. Technical details are provided in Appendix B which contains a paper presented 
by Professor Boyd at the AIAA Joint Propulsion Conference held in Cleveland, July 1998. 
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Detailed Computation of Ultra-Violet Spectra in Rarefied Hypersonic Flow 

Koffi K. Kossif and Iain D. Boyd* 

Cornell University, Ithaca, New York 14853 

A new nonequilibrium radiation model is described that predicts the ultra-violet spectrum of hydroxyl 
(OH) from the electronic transition A—>X for flow conditions corresponding to the Bow-Shock Ultra-Violet- 
2 flight experiment. Unlike previous studies, the new model includes the direct analysis of the electronically 
and vibrationally excited states, OH(A,v’=0,l,2) of hydroxyl in the flow field simulations. The flow field 
is analyzed using the direct simulation Monte Carlo method. Results are presented for the altitude range 
from 80 to 100 km, where the Knudsen number varies from 0.036 to 1.3. The computation uses algorithms 
that improve the numerical resolution of the excited states which typically occur at mole fractions of 10 -15 . 
The collisional transfer of rotational and vibrational energies of the electronic state, OH(A), are also studied 
in detail. It is demonstrated that the usual assumption made in continuum radiation models that the 
rotational and vibrational temperatures of the electronically excited state are the same as those of the 
ground state of the bulk flow fails. An important improvement is achieved in the spectral prediction using 
the new nonequilibrium radiation model and good agreement is obtained between flight data and emission 
predictions over a range of altitudes. 

Nomenclature 

A = excited electronic state designation or reaction rate constant, m 3 molecule -1 s -1 
c = speed of light, m s -1 

E = activation energy or energy level of state, J 

F’ = rotational energy, m -1 

g = statistical weight 

G’ = vibrational energy, m -1 

h = altitude, km 

h = Planck’s constant, 6.625xl0 -34 , J s 

J’ = rotational quantum number 

k = rate coefficient, m 3 molecule -1 s -1 

k = Boltzmann constant, 1.38xl0 -23 , J K -1 

M = a third body species 

n = species number density, m -3 

N = species number of particles or temperature exponent in rate coefficient expression 

f Research Associate. Mechanical and Aerospace Engineering. 

* Associate Professor. Mechanical and Aerospace Engineering. Member AIAA. 
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Q = 

partition function 

t 

time, s 

T t = 

translational temperature, K 

T r = 

rotational temperature, K 

T v = 

vibrational temperature, I< 

V 

cell volume, m 3 

w = 

particle weight 

X 

ground electronic state designation or mole fraction 

e 

vibrational characteristic temperature, K 

v = 

frequency, Hz 

p = 

mass density, kg m~ 3 

r = 

life time or vibrational relaxation time, s 

Subscripts 


e = 

electronic index 

i = 

particle index 

s = 

species index 

r = 

rotation index 

t = 

translation index 

u — 

electronic upper state index 

v = 

vibration index 

V* 

upper vibrational state 

X = 

colliding particle 

0,1,2 = 

excited states index 

oo — 

initial condition 

Superscripts 

diss = 

dissociation 

ex = 

exchange 

form = 

formation 

qu = 

quenching 

N = 

concentration annotation 


Introduction 

Ultraviolet diagnostic experiments have been carried out to measure the intensity and spectral distri¬ 
bution of the radiation emitted from the shock wave of a missile flying at 5.1 km/s (BSUV-2, Ref. 1) and to 
validate aerodynamics and radiation transfer models related to reentry conditions. Detailed spectral mea¬ 
surements from BSUV-2 are available from about 110 km to 71 km with uncertainty of a factor of about 2. 1 
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The Knudsen number range for flow over the 10 cm radius nose of BSUV-2 varies from about 1 to 0.001. 
Therefore, numerical simulation of these flows is difficult due to the thermochemical nonequilibrium condi¬ 
tions and the emission predictions signatures can be in error by orders of magnitude. In a previous article, 2 
we discussed the prediction of ultra-violet emissions radiated by hydroxyl (OH) for hypersonic nonequilib¬ 
rium flow conditions corresponding to BSUV-2. 1 The study obtained flow field solutions which include the 
direct analysis of the electronically excited state species of OH using the direct simulation Monte Carlo 
method (DSMC). Then, the emission was predicted from the flow field solutions using the nonequilibrium 
radiation code NEQAIR. 3,4 The present investigation continues these studies with development of a reliable 
method for predicting the signatures, with the goal of breaking new ground by including the production and 
quenching of the electronically and vibrationally excited states, OH(A,v’=0,l,2) directly in the DSMC flow 
field simulations. The effect of rarefaction, and the extremely small mole fractions of the radiating species, 
about 10“ 9 , have required the development of new algorithms that improve the numerical resolution of the 
trace species. 2 A further goal of the present study is to take advantage of the direct analysis of the electronic 
state species to study in detail its rotational and vibrational energies. This allows assessment of the usual 
assumption made in NEQAIR that the rotational and vibrational temperatures of the excited state species 
are the same as those of the bulk flow or of the ground state species. 

In order to perform the spectral computation, one must first know the number density of the electron¬ 
ically, vibrationally, and rotationally excited states. In the case of an equilibrium condition, this number 
density is determined from the pressure and the enthalpy. However, for nonequilibrium flow conditions, such 
a procedure is not possible. Under some conditions, the number densities of the excited states may instead 
be determined using the Boltzmann distribution for the electronically, vibrationally and rotationally excited 
states. In addition, as described in Ref. 5, the number density of the electronically excited states can also 
be calculated by solving the Master equations. The main assumption in this approach is that a quasi-steady 
state (QSS) exists for the number density of the excited species. The assumption requires the electronic 
excitation rate coefficients to be much faster than the time scales for diffusion and for changes in overall 
properties. Under such conditions, the nonequilibrium flow properties, which are the species number densi¬ 
ties, the rotational, the vibrational and the electronic temperatures, must be known. Then, the vibrationally 
and rotationally excited states number densities are computed using the Boltzmann distribution within these 
levels. However, the present study departs from this approach by computing directly the nonequilibrium 
number density of the electronically and vibrationally excited states species. This approach allows the direct 
assessment of the validity of the QSS theory, especially in high altitude regimes where the chemical processes 
become rare. The present method can also be very useful in helping to understand the nonequilibrium 
processes that the excited species undergo by providing information about their rotational and vibrational 
relaxation phenomena. In fact, in NEQAIR, the vibrational, rotational and electronic temperatures of the 
bulk flow or of the ground state species are used to compute the number density of the electronically, ro¬ 
tationally and vibrationally excited states. In reality, in the conditions of highly nonequilibrium flows, the 
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internal temperatures of the ground state species or those of the bulk flow may be very different from the 
temperatures of the excited states. In such conditions, the approximation used so far in NEQAIR can lead 
to errors by orders of magnitude. In the case of the OH molecule, it was found that the ground state species 
undergoes vibrational energy transfer collisions with rates slower than for the A state, OH(A), by about a 
factor of 100. 6 Therefore, the main goal of the present study is to develop reliable method for signatures 
predictions at high altitude conditions. 

The outline of the present study is as follows: First, the physical models and numerical algorithms 
developed for the DSMC analysis are discussed. Then, the nonequilibrium radiation code used in this study 
is briefly described. The model of vibrational energy transfer of the electronically excited state is introduced. 
In terms of results, general properties of the DSMC flow field solutions are presented. Results obtained with 
the new nonequilibrium radiation model are compared with the NEQAIR model. Direct comparison of the 
predicted peak emissions and spectra with the BSUV-2 data are presented. Finally, the flow field properties 
that include the rotational and vibrational temperatures of OH(A) are presented and the emission predictions 
are analyzed in terms of the OH(A) rotational temperature dependence. 

Numerical models 

The numerical models employed in the present work are described in this section. The physical models 
implemented in the DSMC code are first considered. The new nonequilibrium radiation model is then 
introduced and finally, the vibrational energy transfer of OH(A) is presented. 


Physical models 

The DSMC code employed in the present study is based on the vectorized algorithm described in 
Ref. 7. The code includes finite rates of rotational and vibrational relaxation, and dissociation, exchange, 
and recombination chemical reactions. A key development in the code is the implementation of a numerical 
weighting scheme that allows accurate and efficient simulation of trace species. 8 For the altitude range 
of interest here, the weighting scheme is needed to resolve the excited states, OH(A,v =0,1,2) and their 
precursors, OH and H 2 0. It will be found from computations that the electronically and vibrationally excited 
state species, OH(A,v’=0,l,2) exist in the flows at mole fractions of less than 10~ 9 and these conditions are 
handled efficiently by the new weighting scheme. However, in order for any OH(A,v’=0,l,2) particles to 
exist in the flows, they must be produced through collisional excitation. For the altitudes of interest here (80 
to 100 km), there are few intermolecular collisions and therefore the occurrence of possible collision events 
is very rare. As described in Ref. 2, whenever a collision involving an OH particle has sufficient energy to 
allow electronic and vibrational excitation, the reaction probability is set to unity. The numerical weights 
of the products are then adjusted to ensure conservation of mass. A special procedure is employed for 
these excitation collisions whereby several OH(A) particles are created instead of just one. Once again, the 
numerical weight of each particle is adjusted to conserve total mass. Using this procedure, it is possible to 
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maintain thousands of OH(A,v’=0,l,2) particles in the simulation. The spontaneous de-excitation reactions: 

OH(A,v' = 0,1,2)^ OH + hv (1) 

and the predissociation reaction: 

OH(A,v' = 2)^0 + H (2) 

are implemented in the DSMC code using the algorithm described in Ref. 2. 

The observed OH ultra-violet emission is from the A—>X transition only. Therefore, it is relatively 
simple to include the production and quenching of the electronically and vibrationally excited states directly 
in the DSMC flow field simulations. An 11-species reacting flow model is employed for the following species: 
N 2 , N, O 2 , O, NO, and H 2 O, OH, H, and electronically and vibrationally excited states OH(A,v’=0,l,2). 
The maximum vibrational levels in the computations is fixed to v’=2 because the experimental spectra 
show no features of vibrational levels greater than 2. In addition to the air reactions, 7 the hydrogenated 
species undergo the reactions listed in Table 1. Note that the quenching and de-excitation coefficient rates 
are from Ref. 4 and the reverse reactions coefficient rates are then computed using the principle of balance 
at equilibrium. The reactions are simulated by the DSMC code using the total collision energy form of 
the Generalized Collision Energy model. 9 In previous studies, 2 ' 10 it was shown that prediction of ultra¬ 
violet emissions was sensitive to the model assumed for interaction of the gas with the vehicle surface. 
It was found in Ref. 10 that the best agreement of prediction with experimental data occurred for the 
translational accommodation coefficient of 0.85. This value is consistent with molecular beam experiments. 
The rotational and vibrational accommodation coefficients were then set respectively to 0.5 and 0.1. The 
rotational accommodation coefficient is also consistent with observed molecular beam, and the vibrational 
accommodation coefficient was estimated. In the present work, these values are retained and the surface 
temperature is taken as 500 K as measured during the BSUV-2 flight. It is also shown in Ref. 2 that the 
prediction of ultra-violet emission for OH was sensitive to the free stream concentration. Therefore, the 
concentrations of the hydrogenated species obtained from a model atmosphere, 11 and used in Ref. 2 are 
retained here and listed in Table 2. 

Nonequilibrium Radiation model 

The nonequilibrium radiation code described in Ref. 3 and employed for prediction of emission from flow 
field solutions, was discussed and verified in Ref. 4 for OH comparing the measured spectra with the predicted 
emissions. The common assumption made in the NEQAIR code is that the quasi-steady-state (QSS) exists for 
the number densities of the electronically excited species and a single rotational and vibrational temperature 
model is employed. The validity of the QSS assumption for high altitude flows was assessed in Ref. 2 using 
the number density of electronic state species from the DSMC computation. In that approach, the number 
density of the rotationally and vibrationally excited state was computed using the Boltzmann distribution. 
The number density of the upper state in the spontaneous emission modeling is then written as 
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(3) 
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where n e represents the number density of the electronic state. The present work, however, departs from 
this approach by calculating directly from the DSMC code the number densities of the electronically and 
vibrationally excited states. In this new approach, the number density of the upper state is written as 


(2J 1 +1) f hc.F'(J ') yj 


(4) 


In addition to assessment of the QSS assumption, the new model allows direct assessment of the usual 
assumption made in NEQAIR that the rotational and vibrational temperatures of the excited state are the 
same as those of the bulk flow or of the ground state species. 


Rotational and vibrational energy transfer of OH(A) 

Collisional transfer of rotational and vibrational energies of OH(A) are the topic of interest in this 
section. Simulating directly in the DSMC code the production and quenching of OH(A), provides the 
opportunity to study in detail the rotational and vibrational energy transfer. Prom the exchange reaction: 


OH (A, v' = 1) + M -* OH(A, v‘ = 0 ) + M (5) 

where, M = Ar 2 ,02 

we can define the vibrational relaxation time, t as 

1/r = n x Kio(l - exp{-O v /T) ~ n x Kio0 v /T ( 6 ) 

where K w is the vibrational excitation rate related to the v’ = 1 -» 0 vibrational transition cross sections 
which are reported in Ref. 4 and listed in Table 1. The vibrational relaxation times of OH(A) and OH with 
N 2 as a colliding partner are shown in Fig.l. It is noted that the vibrational relaxation time of OH(A) is 
orders of magnitude lower than that of OH which indicates the OH vibrational energy transfer rate to be 
slower than for the A state, OH(A). The rotational relaxation times of OH and OH(A) are assumed to be 
the same. 

The rotational and vibrational energy transfer are then implemented in the DSMC code in addition to the 
other molecules using the algorithms described in Refs. 12 and 13. However, for vibrational energy transfer 
modeling, care must be taken when simulating both vibrational energy exchange and chemical exchange 
reactions between OH(A,v’=0,l,2). In fact, the exchange reaction processes involving OH(A,v’=0,l,2) also 
characterize the vibrational energy exchange phenomena. Since the vibrational levels of OH(A) are each 
treated as if they were separate chemical species and the vibrational relaxation is then treated as a chemical 
reaction, no additional vibrational relaxation process is applied to OH(A). 

Results 
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The presentation of results is divided into three sections: (1) variation of flow field properties over the 
altitude range of interest; (2) comparisons of predicted emissions using the new nonequilibrium radiation 
model with the NEQAIR model predictions and BSUV-2 flight data; and (3) presentation of general flow 
properties including OH(A) rotational and vibrational temperatures and sensitivity of predicted emissions 
to OH(A) rotational temperature. 

As in the previous computational studies of the BSUV-2 flight, only flow over the 10 cm radius spherical 
nose of the vehicle is simulated. The velocity is in every case 5.1 km/s, and four different altitudes are 
considered: 80, 88, 94, and 100 km. The hydrogenated species mole fractions, obtained from a model 
atmosphere presented in Ref. 11. are listed in Table 2. In all cases, the accommodation coefficients for the 
three energy modes (translation, rotation and vibration) are set respectively to Ot—0.85, ov=0.5, and q v = 0.1 
and the surface temperature is taken as 500 K as measured during the BSUV-2 flight. It is also necessary 
to note that the values of the accommodation coefficients of the internal energy modes do not influence the 
radiation prediction. 

Flow field Properties 

The computed temperatures and number densities are shown in Figs. 2a through 3b for altitudes of 80 
and 100 km respectively. The results show the shock layer becoming thicker at higher altitudes and that 
a large degree of thermal nonequilibrium exists for the altitudes investigated. At 80 km, the extent of the 
chemistry for some of the specie.^ is illustrated in Fig. 2a where OH and OH(A,v'=0,l,2) number densities 
profiles are included. However, at 100 km altitude, the number densities of species are almost flat indicating 
the relative importance of diffusion over chemical processes and the peak number density of OH(A,v’=0) is 
reduced from 10 12 m -3 at 80 km to about 10 7 m -3 at 100 1cm. The peak values of the temperatures are 
also considerably reduced indicating the effect of rarefaction. The effect of rarefaction is also illustrated in 
Fig. 3a by the OH(A,v’=2) number density profile. In fact, for any OH(A,v’=2) particle to exist in the flow, 
it must be produced through OH( A,v’=0) or OH(A,v’=l) collisional excitation. Despite the fact that several 
OH(A,v’=2) particles are assumed to be created at each collisional event, very few OH(A,v’=2) particles 
are simulated after several thousand time steps. Note that, despite the efficiency of the numerical model in 
handling the trace species simulation, the statistical noise related to OH(A) properties are significant. For 
clarity of presentation, the shock layer domain is reduced in Figs. 3a and 3b. This explains the fact that the 
temperatures do not reach the free stream value. 

In Figs. 2b and 3b, the vibrational temperature of OH(A) is computed from the predicted number 
density, n Q of the ground state, OH(A,v’=0), and the first excited state, OH(A,v’=l) number density, n\ 
using the Boltzmann distribution: 


T„ = — 


AEJk 

"■(-) 
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where A E v characterizes the difference between the two energy levels. It is noted that the vibrational 
temperature computed in this manner is relatively constant. 

OH Emission 

In Figs. 4a and 4b, the computed spectra at 80 and 100 km are compared with the results obtained 
using the QSS theory. It is shown that the spectra computed using the QSS theory do not simulate the 
vibrational transition of OH(A) 1-0 at 0.28 /xm and the shoulder of the OH(A) 0-0 vibrational transition at 
0.31 /xm is too narrow. The absence of the 1-0 feature in the QSS theory is explained by the inaccuracy of the 
vibrational temperature used in this model. Note that the new radiation model does not use the vibrational 
temperatures, rather the number densities computed from the flow field simulations. The splitting of the 0-0 
peak and the left shoulder of the 0.31 /im peak depend on the rotational temperature, and the disappearance 
of the 0-0 peak at 100 km in both models is due to the small value of the rotational temperature. 

In Fig. 5, the spectrum computed with the new model is compared with the BSUV-2 flight data at 
100 km. While the comparison is reasonable, the peak value at 0.28 /xm in the prediction is relatively small 
compared to the data and the left shoulder of the 0.31 /xm peak has almost disappeared. One possible 
explanation for the relatively small value of the predicted peak at 0.28 /xm is inaccuracy in the OH(A) 
collision cross sections. In fact, the ratio of the peak heights at 0.28 /xm (v' = 1 —* v' = 0) and 0.31 /xm 
(v' — 0 -+ v' = 0) depends on the vibrational temperature. The collision cross sections for vibrational 
energy exchange reactions, listed in Table 1, were measured for temperatures of about 3,000 K and the 
use of these data for flow conditions with translational temperatures around 10,000 K may lead to errors. 
The disappearance of the left shoulder of the 0.31 /xm peak is explained by the relatively small value of the 
computed rotational temperature indicating that the simulation cannot predict correctly all spectral features 
of the data. As shown in Fig. 3b, the rotational temperature is about 1,500 K, whereas the spectral data 
indicate the rotational and vibrational temperatures to be on the order of 3,000 K and 7,000 K respectively. 

The peak emissions at 0.31 /an as a function of altitude are compared in Fig. 6. The experimental data 
are from the spectra measured by BSUV-2. The spectral peak height at 80 km is derived from calculation 
and measurement data using physical arguments. 4 The agreement is quite good at all altitudes lying within 
a factor of 2 at 88 km and about 10 at 100 km. Recall that there is about a factor of 2 uncertainty in the 
BSUV-2 data. 

OH(A) Flow Properties and Emission Sensitivity to Rotational Temperature. 

To assess the validity of the single rotational and vibrational temperatures assumption employed in 
NEQAIR, profiles of bulk flow and OH(A) temperatures along the stagnation streamline are shown in 
Figs. 7a and 7b for altitudes of 80 and 100 km. Note that the vibrational temperatures obtained through 
vibrational relaxation of OH(A) and reproduced in Figs. 7a and 7b, are different from the temperatures 
in Figs. 2b and 3b, which are obtained assuming a Boltzmann distribution. It is seen in Fig. 7b that the 
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internal temperatures of OH(A) are very high in the free steam. This is explained by the importance of the 
diffusion process. In fact, due to rarefaction effects at 100 km, OH(A) particles produced through chemical 
processes diffuse all along the stagnation streamline from their point of formation preventing OH(A) internal 
temperatures from converging to the free stream value. It is also clearly shown that the internal temperatures 
of OH(A) are different from those of the bulk gas. However, at 80 km where there is a sufficient number 
of collision events, the two rotational temperatures are very similar. The two vibrational temperatures 
are quite different; the OH(A) vibrational temperature is very close to its rotational temperature and the 
peak value is about 7,500 K. At 100 km, simulations show that OH(A) internal temperatures are very high 
compared to the bulk gas. The OH(A) rotational and vibrational temperature peaks are about 4,000 K and 
5,000 K, whereas the bulk gas flow simulation indicates the flow rotational and vibrational temperatures 
to be about 1,500 K and 800 K respectively (Fig. 3b). The present results also show that the OH(A) 
vibrational temperature becomes important compared to the rotational temperature in agreement with the 
BSUV-2 flight experiments. However, the present results are different by about 1,000 K from the spectral 
data which indicate the OH(A) rotational and vibrational temperatures to be on the order of 3,000 K and 
7,000 K respectively. To explain the discrepancy, further results were computed by varying the internal 
energy mode accommodation coefficients. However, the results obtained for several tests indicate that the 
internal temperatures do not vary much with accommodation coefficients. It is concluded that the internal 
temperatures of OH(A) are related to OH(A) chemical processes and this conclusion will be confirmed in the 
following discussions. In simulating OH(A) production from OH molecules using the generalized collision 
energy model, any excess collision energy above the activation energy is equipartitioned into the rotational 
and vibrational modes. This approach is partially verified through the relatively good agreement obtained 
with the BSUV-2 spectra. 

Sensitivity of the emission predictions to the OH(A) rotational temperature is considered here. The 
new nonequilibrium radiation model is again considered. In contrast to results shown above, the rotational 
temperature of the bulk gas in the emission prediction is replaced by the OH(A) rotational temperature. 
As seen from Fig. 8b at 100 km, the splitting of the 0-0 peak and the right shoulder of 0.31 fim peak vary 
considerably with the rotational temperature. At 80 km (Fig. 8a), this variation is insignificant indicating 
that the OH(A) rotational temperature is the same as the bulk rotational temperature. In Fig. 9a, the 
spectrum computed using the OH(A) rotational temperature is compared with the BSUV-2 data at 100 km. 
An improvement is obtained with the left shoulder of the 0.31 ^rm peak which has increased significantly 
in the calculated spectrum. However, this peak value is relatively high compared to the experimental data, 
which indicates that the OH(A) rotational temperature is too high in the numerical simulation. To assess 
this conclusion, the emission is predicted at 100 km altitude assuming a rotational temperature of 3,000 K. 
Figure 9b indicates a decrease of the right shoulder of 0.31 gm peak which is now in good agreement with the 
BSUV-2 flight data. This result confirms the OH(A) rotational temperature to be on the order of 3,000 K. 
Using the OH(A) rotational temperature, it is noted that the peak emission at 0.31 /im is reduced, especially 
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at 100 km where the OH(A) ancl the bulk gas rotational temperatures have a large discrepancy. 

Assuming the rotational temperature of OH(A) to be 3,000 K at 100 km altitude and comparing the 
predicted peak with the flight data, an estimate of the OH free stream mole fraction is then obtained. The 
derived OH mole fraction is about 4.0xl0 -9 . This value is about 20 times higher than the value reported 
in Table 2. Usiiig the same procedure, the free stream water mole fraction is derived at 80 km through 
comparison of the predicted peak with the derived 80 km data point. The derived water mole fraction is 
about 6.1xl0~ 6 and very close to the data reported in Table 2. It was found in Ref. 2 that at low altitudes 
such as 80 and 88 km, the emission predictions were directly proportional to the ambient water mole fraction. 
However, at high altitudes, such as 94 and 100 km where the effects of rarefaction are very significant, the 
emission predictions were directly proportional to the OH free stream mole fraction. Considering these 
results, the ambient water mole fraction and the OH free stream mole fraction are also derived at 88 and 
94 km. The water mole fraction is found to be on order of 4.2xl0 -6 at 88 km and the OH mole fraction is 
about 6.0xl0 -9 at 94 km. The derived mole fractions of the hydrogenated species are reported in Table 3. 

Using the mole fractions of the hydrogenated species listed in Table 3, the DSMC method is applied 
to compute the flow field properties, and the new nonequilibrium radiation model is used to predict OH 
emission over the altitude range of 80 to 100 km. In h ig. 10, the peak emissions at 0.31 /rm as a function of 
altitude are compared with data taken from the spectra measured on the BSUV-2 flight. Clearly, excellent 
agreement is obtained suggesting that the hydrogenated species mole fractions experienced in the BSUV-2 
flight experiments are within the data listed in Table 3. However, as mentioned previously, the quenching 
coefficient rates can also explain the discrepancy observed between the predicted spectra and the experimental 
data. The derived mole fractions must then be considered representative rather than definitive. 

Considering once again the flow and OH(A) properties at 100 km altitude, it is shown in Fig. 11 that 
the OH(A) rotational and vibrational temperatures are in good agreement with the temperatures indicated 
by the BSUV-2 spectral data. It is noted that the internal temperatures of OH(A) obtained with the new 
free stream mole fraction are higher than those reported in Fig. 7b. In fact, since the OH mole fraction is 
increased, the occurrence of collisional excitation is also increased and so are the number of particles in high 
rotational and vibrational levels. This result also demonstrates that, in high altitude conditions, OH(A) 
internal energy temperatures are related to the chemical processes rather than collisional energy transfer. 

Summary 

The direct simulation Monte Carlo method was used to predict the number densities of electronically 
and vibrationally excited state species of OH in nonequilibrium, hypersonic flight conditions. The conditions 
corresponded to the Bow-Shock Ultra-Violet-2 flight experiment over the altitude range of 80 to 100 km. 

The profiles of OH(A,v’=0,l,2) were then used lo predict the ultra-violet spectra from OH(A). The 
comparisons of both the peak value and the overall spectral shape were found to be in good agreement with 
the flight data. However, the ratio of the peaks at 0.28 jum showed some disagreements which could be 
related to the inaccuracy of the quenching and exchange collision cross sections at high temperature. In 
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addition, the left shoulder height at 0.31 pm showed disagreements at high altitude which are due partly to 
inaccuracy in computing rotational temperatures. Nevertheless, the present study represents a significant 
improvement compared to NEQAIR. 

Comparisons of OH(A) rotational and vibrational temperatures with the bulk gas flow values showed 
that at the lower altitude conditions, the two rotational temperatures were in good agreement while the 
vibrational temperatures were very different. Furthermore, it was also shown that at high altitude, the 
OH(A) internal temperatures were higher than the bulk gas. Using OH(A) rotational temperatures, the 
emission was computed that reproduced both the spectral shapes and the peak emissions dependence on 
the rotational temperature. Comparisons of the emission predictions with the BSUV-2 flight experiments 
allowed the estimation of H 2 0 and OH ambient mole fractions over the range of altitudes of interest. The 
H 2 0 ambient mole fractions were found to be close to values from a model atmosphere, while OH free stream 
mole fractions were found to be higher by about a factor of 20. Using the derived free stream mole fractions, 
the flow field was then computed and the emissions predicted. Good agreement was obtained with the flight 
data for both the peak emission and the overall spectra. With the new free stream mole fractions, agreement 
was obtained between OH(A) internal temperatures and those indicated by the spectral data. 
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Figure Captions 


Fig. 1 OH and OH(A) vibrational relaxation times as a function of temperature. 

Fig. 2: a) Number densities and b) temperatures along the stagnation streamline at 80 km altitude. 

Fig. 3: a) Number densities and b) temperatures along the stagnation streamline at 100 km altitude. 

Fig. 4 Comparison between DSMC and NEQAIR of the ultra-violet spectra: a) at 80 km altitude and b) at 
100 km altitude. 

Fig. 5 Normalized ultra-violet spectra at 100 km altitude. 

Fig. 6 Peak emission at 0.31 fim as a function of altitude. 

Fig. 7 Temperatures along the stagnation streamline: a) at 80 km altitude and b) at 100 km altitude. 

Fig. 8 Rotational temperature dependence of the ultra-violet spectra: a) at 80 km altitude and b) at 100 km 
altitude. 

Fig. 9 Comparison of the normalized ultra-violet spectra at 100 km: a) OH(A) rotational temperature and 
b) rotational temperature of 3000 K. 

Fig. 10 Peak emission at 0.31 fin i as a function of altitude using the derived OH mole fractions. 

Fig. 11 Temperatures along the stagnation streamline at 100 km altitude using the derived OH mole fractions. 
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Table 1 OH kinetics mechanism 


Reaction 

Rate coefficient, 
m 3 molecule"* 1 s" 
k=AT N e(- E ^ 

A 

-1 

N 

E/Jfc.K 

H 2 0 + N 2 ->• OH + H + N 2 

5.81xl0 -15 

0.00 

-53,000.0 

H 2 0 + 0 2 -*■ OH + H + 0 2 

1.13xl0-° 7 

-1.31 

-59,400.0 

h 2 o + o-+oh + h + o 

1.13xlO -07 

-1.31 

-59,400.0 

OH + N 2 -» 0 + H + N 2 

1.25xl0 -15 

0.06 

-51,000.0 

OH 4- O 2 —* 0 + H 4- O 2 

1.25xl0 -15 

0.06 

-51,000.0 

0H + 0-»0 + H + 0 

1.25xl0~ 15 

0.06 

-51,000.0 

H + 0 2 4 —♦ OH + 0 

3.65xl0 -16 

0.00 

-8,450.0 

0 + H 2 0 <—- OH + OH 

1.13xl0 -16 

0.00 

-9,240.0 

OH + N 2 -* OH(A,v’=0) + N 2 

9.40xl0 -17 

-0.37 

-46,600.0 

OH + 0 2 -> OH(A,v’=0) + 0 2 

2.70xl0 -15 

-0.37 

-46,600.0 

OH + N 2 —> OH(A,v’=l) + N 2 

2.30xl0 -19 

0.27 

-51,173.0 

OH + 0 2 —> OH(A,v’=l) + 0 2 

6.68xl0~ 18 

0.27 

-51,173.0 

OH(A,v’=0) + N 2 - OH + N 2 

1.64xl0 -19 

0.50 

687.0 

OH(A,v’=0) + 0 2 —> OH 4- 0 2 

4.76xl0~ 18 

0.50 

450.0 

OH(A,v’=l) + N 2 — OH + N 2 

2.17xl0~ 19 

0.50 

687.0 

OH(A,v’=l) + 0 2 —> OH + 0 2 

6.31xl0 -18 

0.50 

450.0 

OH(A,v’=0) + N 2 -» OH(A,v’=l) + N 2 

5.65xl0~ 18 

0.50 

-4,573.0 

OH(A,v’=0) + 0 2 - OH(A,v’=l) + 0 2 

9.66xl0“ 19 

0.50 

-4,573.0 

OH(A,v’=0) + N 2 OH(A,v’=2) + N 2 

7.46xl0 -18 

0.50 

-9,146.0 

OH(A,v’=0) + 0 2 -» OH(A,v’=2) + 0 2 

2.43xl0 -19 

0.50 

-9,146.0 

OH(A,v’=l) + N 2 - OH(A,v’=2) + N 2 

S.OlxlO -18 

0.50 

-4,573.0 

OH(A,v’=l) + 0 2 -♦ OH(A,v’=2) + 0 2 

9.20xl0 -19 

0.50 

-4,573.0 

OH(A,v’=l) + N 2 -> OH(A,v’=0) + N 2 

5.37xl0~ 18 

0.50 

0.0 

OH(A,v’=l) + 0 2 -» OH(A,v’=0) + 0 2 

9.17xl0~ 19 

0.50 

0.0 

OH(A,v’=2) + N 2 -4 OH(A,v’=0) + N 2 

6.71xl0 -18 

0.50 

0.0 

OH(A,v’=2) + 0 2 -» OH(A,v’=0) + 0 2 

2.18xlD -19 

0.50 

0.0 

OH(A,v’=2) + N 2 -» OH(A,v’=l) + N 2 

7.61xl0 -18 

0.50 

0.0 

OH(A,v’=2) + 0 2 -» OH(A,v’=l) + 0 2 

8.74xl0~ 19 

0.50 

0.0 

Reaction 

Einstein coefficient, s 1 


OH(A,v’=0) -» OH + hv 

1.44x10® 



OH(A,v’=l) —> OH + hi/ 

1.36x10® 



OH(A,v’=2) -+ OH + hi/ 

1.35x10® 



Predissociation reaction 

rate, s"" 1 



OH(A,v’=2) -> 0 + H 

3.70x10® 
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Table 2 Flow conditions 


h (km) 

Poo (kg m 3 ) 

Too (K) 

X// 2 o 

XoH 

Xh 

80 

2.00xl0 -5 

181 

5.6x10"® 

4.3xl0 -9 

2.0xl0" 7 

88 

5.11xl0~® 

195 

3.0x10-® 

1.2X10" 9 

7.0x10-® 

94 

1.30x10-® 

177 

1.2x10-® 

3.0xl0 _1 ° 

9.0x10"® 

100 

5.68xl0~ 7 

185 

7.2xl0 -7 

2.0xl0 _10 

l.OxlQ -5 


Table 3 Derived freestream mole fractions. 


h (km) 

Xh 2 o 

XoH 

80 

6.1x10-® 


88 

4.2x10-® 


94 


6.0xl0- 9 

100 


4.0xl0 -9 
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Abstract 

The direct simulation Monte Carlo method is employed to compute cold flows of ammonia for a large 
arcjet that is to be tested in space in the upcoming ESEX flight experiment. The nozzle flow computation 
indicates that the flow is almost in thermal equilibrium at the nozzle exit. A very large and expensive com¬ 
putation of the back flow region of the actual spacecraft geometry is performed to provide predictions of mass 
fluxes that will be measured in flight by quartz crystal micro-balances. It is indicated that contamination 
of the spacecraft occurs even in regions lying behind a plume shield. A further computation is performed 
to simulate the interaction of the arcjet plume with the ambient atmosphere. The high impact energy is 
offset by the very low atmospheric density at the spacecraft operational altitude of 833 km. Nevertheless, it 
is indicated that ammonia chemistry occurs and the primary products are NH, NH 2 , and OH. These species 
radiate strongly in the ultra-violet. The estimated emission intensities of the molecules are similar to those 
measured previously in situ by a lower velocity reentry experiment. An estimate is also made of the intensity 
of emission from ammonia. In all cases, it is concluded for the cold flow that none of the emissions will be 
detectable by the ground based observation facility that is part of the space experiment. 

Introduction 

Arcjets represent a mature electric propulsion technology that is replacing chemical propulsion engines 
for orbit maintenance procedures on many satellites. The U. S. Air Force has developed the Electric Propul¬ 
sion Space Experiment (ESEX) 1 to study several of the effects of using a very high power ammonia arcjet on 
a spacecraft. The ESEX experiment is part of the ARGOS flight that will orbit the Earth at an altitude of 
about 833 km. The flight experiment, due for launch in the winter of 1998, includes diagnostic measurements 
to analyze exhaust plume contamination and radiative heating effects as well as arcjet performance. Remote 
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3 Post-doctoral Research Associate. Mechanical and Aerospace Engineering. 

4 Research Professor. Chemistry. Member AIAA 
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optical measurements are also planned from a ground station at Maui. Several arcjet operational modes will 
be investigated including the cold flow case where the arc is not ignited. 

The goal of this paper is to conduct numerical simulations of the cold flow of the arcjet. The primary 
motivation for these computations is to provide estimates of the contamination on the spacecraft by ammonia 
that will be measured by quartz crystal micro-balances. Hence, the focus of these computations will be the 
back flow region behind the thruster exit plane. A secondary focus of the study is to investigate the interaction 
of the ambient atmosphere with the ammonia jet exhausting from the thruster. In particular, the goal is to 
estimate the radiative emission generated under these conditions. 

The layout of the paper is as follows. First, the arcjet thruster and diagnostics of ESEX are described. 
Next, the computational approach is discussed. Due to the relatively low densities involved, the direct 
simulation Monte Carlo method (DSMC) 2 is employed. Discussion of results is divided into separate sections 
concerning the back flow contamination studies and the radiative emission calculations. 

The ESEX Experiment 

A schematic diagram of the ESEX hardware is shown in Fig. 1. As described in Ref. 1, ESEX consists 
of a 26 kW ammonia arcjet with associated hardware and a diagnostics package. The arcjet has a throat 
radius of 1.9 mm and an exit radius of 12.6 mm. The nozzle half angle is 19°. The nominal flow rate is 
0.25 g/s of ammonia. At full power this provides a specific impulse of 800 s and a total thrust of 2 N. 

The ESEX hardware, including the diagnostics package, is part of the ARGOS spacecraft. The diagnos¬ 
tics will measure several aspects of thruster operation including: (1) performance, from thrust measured by 
an accelerometer, (2) thermal radiation, measured using radiometers, and (3) electromagnetic interference 
using two spiral cavity antennas. Of most interest to the present study is the measurement of exhaust plume 
contamination using quartz crystal micro-balances (QCM’s). These instruments provide species-resolved 
mass flux data. The QCM’s are located in four different positions. QCM’s 2 and 3 are located on the diag¬ 
nostic platform in the back flow region and their inlet planes are parallel to the thruster exit plane. QCM 1 
is located on a solar cell witness plate and points towards the thruster. QCM 4 is located on a deployable 
boom that extends beyond the perimeter of the spacecraft and again its inlet plane is parallel to the thruster 
exit plane. 

The emissions radiated by the arcjet during ignited operation will be measured from a ground-based 
facility at Maui. From laboratory experiments, the principal spectral features easily observable from the 
ignited arcjet are expected to arise from electronically excited NH (in the 1st excited, or A, state) in the 
near-UV at 330-345 nm, and from the Balmer series of hydrogen in the red, blue, and near-UV. 3 

Numerical Simulations 

The mass flow rate and physical dimensions of the ESEX thruster indicate that the nozzle flow is in the 
near continuum regime with a Knudsen number at the throat of about 10 4 . By the nozzle exit plane, the 
Knudsen number has increased to about 10~ 2 . In most of the plume expansion region the Knudsen number 
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« , ,v : s dearly in the rarefied regime. Therefore, all of the computations 

te si6ni8c ““ y h ' 6 eI “ stalato Monte Carlo method (DSMC). The DSMC technique 

Nozzle Flow * 

The nozzle Sow is computed separately from the plume expansions to reduce computational exp • 

DSMC nozzle code from previous resisto-jet studies' is employed again here. The computation egms J 

le throat and the initial starting condition is obtained from isentropic t eory. 
downstream of the nozzle throat and the. ^ ^ ^ ^ a stagnaliM1 temperature 

necessary properties are obtained from , 1 Q00 b 10 q cells 

of 295 K The computation is performed on a desh-top personal computer using a gnd o, 1.000 10 c 

° , , , . „„ ooo particles The exit plane profiles of number density, temperature, and velocity 

L« in Figs. 2a.2c. They indicate a relatively large boundary layer and thermal —,el 

•n *orm data. obtained from isentropic theory are also included. 

spedficTnqpulse^is 07 ZZ the' total thrust is 0.24 N. These can be compared with the ignited arc*, 
performance of 800 sec and 2.0 N respectively. 

Plume Back Flow 

Accurate plume bach flow computations can only be performed with the DSMC Unique due to he 

•tv In this regime continuum methods simply breakdown, and free molecular, 
rapid decrease in density. In this regime, cun Obtaining 

. from shadowing effects in the complex geometries of real spacecra . 

of-sight techniques suffer from shadowing e . , t tn ffet anv sienif- 

, DSMC tesults for back flow configurations is not easy however since it is difficult ge 
accurate DSMC resu t > stoJ ies, a very large DSMC computation 

ic“, —' P ” ,,dK the “ Zx diagnostic platform is performed using a par-, numerically 
of the arc*, plume expandingP ^ ^ „ Heally soiled 

efficient implementation o tUVSMC ^ ^ ^ ^ ^ ^ tke variation 

for this application since d u. ^ mate # possible „ resolve „ a numerically efficient 

of time steps and particle weigni , ,. , .. thruster The 

manner the large gradients in density experienced as the gas flows into the 

computation begins a. the nozzle exit plane using output from the nozzle compuU ron. Tie 

to make sure the cell size is proportional to the loc* mean free path everywhere. I. ■ -—* 

. . nerfect va cuum The total atmospheric number density at an altitud 
plume expands into a perfect vacuu . -p; o a 

about 10- m-t. This is ten orders of magnitude less than the nozzle exit density as shown Frg. . 
Plume- Atmosphere Interaction 

The ESEX arcjet is configured as to fire directly along the velocity rector into the oncoming^am 
flow. At an altitude of 833 km, the main constituents in the atmosphere a. night are *mx: y ™ ^ 

helium (39%) and atomic oxygen (1%). At orbital speed, the collision energy of an ammo 
either helium or oxygen is sufficient to dissociate the ammonia with the heavier oxygen being significantly 
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more reactive. There are a number of possible chemical channels and some lead to production of NH and OH 
which are both strong radiators in the ultra-violet. All other species except ammonia will be present only 
in trace amounts in the plume flow. Procedures for simulating trace species in the DSMC technique have 
been developed 13 and applied to model OH ultra-violet emission measured in a reentry flight experiment. 
Those procedures are again used here along with the chemical reaction mechanism shown in Table 1. The 
table also lists the sources of the reaction rates. 

The flow from the arcjet begins at the thruster exit plane using the data from the nozzle simulation. 
The ambient atmospheric species of H, He and 0 are introduced at the downstream boundary of the flow 
domain at orbital speed. 

Results 

The only results of interest from the nozzle computation are those at the nozzle exit plane and these 
were presented earlier. As discussed previously, these results are used to start both the contamination and 
plume emissions studies that are now discussed in detail. 


Contamination 

The MONACO computation of the impingement of the ammonia plume onto the ESEX payload is 
performed on 32 processors of an IBM SP-2 computer and uses 11 million particles with 34,000 cells. The 
large size of the simulation is required to ensure that all of the tenuous back flow region is populated wtth 
particles. The grid employed for the computation is shown in Fig. 3a. The locations and orientatrons of the 

four QCM’s are indicated. 

Contours of number density are shown in Fig. 3b. They reveal that the number density decreases by 
7 orders of magnitude as the gas flows from the thruster exit to the diagnostic platform. The contours 
of translational temperature shown in Fig. 3c indicate the regions of expansion and recompression in the 

complex spacecraft geometry. 

Spacecraft surface impingement properties are shown in Figs. 4a and 4b as a function of distance along 
the spacecraft starting at the lip of the nozzle in the exit plane. In Fig. 4a, it is shown that the mass flux 
to the ESEX diagnostics platform decreases by 6 ordem of magnitude. The locations of QCM’s 2 and 3 are 
indicated. Table 2 lists the mass fluxes predicted at all four QCM locations. QCM locations 2 and 3 are on 
the platform surface and the predicted mass fluxes are evaluated directly from particles striking that surface. 
QCM locations 1 and 4 are not located on any solid surface in the simulation and so are evaluated from the 
product of density and velocity resolved in the direction normal to the QCM inlet face. 

In Kef. 1, estimates of the mass fluxes expected at the QCM locations for the ignited arcjet were 
provided. These are also listed in Table 2 for comparison. It is expected that the mass flux for the cold 
flow will be greater in the back flow region than that for the ignited flow. This is because the energy from 
the arc generates very high velocity flow at the nozzle exit creating a well defined, forward movmg beam. 
This conjecture is supported by the predictions for the mass fluxes at QCM locations 1 and 4. For the CO 
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— flow, the mass flux a, QCM 4 is two orders o. magnitude higher than that es.taa.ed incite 

the ignited flow. Surprisingly, the predictions for total mas, flux onto the dragnostrcs plat om (Q 

nTlH o\ f or the cold and ignited flows are almost identical. 

I„ Fig 4b, the surface pressure is show, to decrease by seven orders of magnitude from the nose e rp 

to the diagnostic platform. Clearly, the use of the arc*, in cold flow mode will no. produce any srgnrfican 
momentum transfer to the vehicle. 

Plume Emission 

*. total physical domain of 50 cm aerially and 25 cm radiaUy is considered. A total of 40,000 cells wrth 

500 000 particles is employed. The simulation is performed on a personal computer. 

' contours of the number density of ammonia is shown in Fig. 5a. There is no appreciable srgn tha 
C u . bv the ambient atmosphere. In contrast, the ambient atoms are clearly interacting 

!i7.he“ uld’lr plume. This is also made dear in Fig. 5b tha. shows the translarional and rotottaml 
temperatures of ammonia along the plume axis. Note that the vibrational mode is included in 

but is unexcited throughout the nossle and plume flow. „ 

When the properties of the trace species are investigated, a very different picture of the flow e g ^ 
The level of chemical activity in this highly energetic and rarefied flow is shown in Fig- a w ic P 
some of the species number densities along the plume arris. The He proflle is almost 

and the NH, profile is almost identical to that for OH and so for clarity these ata are no ^ 

the ambient atoms dissociate a small fraction of ammonia molecules into NH„ NH and OH. 
molecules are strong razors in the ultra-violet, so it is appropriate to estimate the emission that m 

ge “Tta 771777Tternal modes of OH are shown along the plume arris in Fig. 6b. The very 
energetic dissociation collisions produce OH molecules that are highly excited. Through ^ 

the translational and rotational temperatures are reduced. However, the vibmtional ~, t 
sufficiently slow tha. almost no further change occurs to the vibrational energy thus earing 

Thp monerties of OH can be compared with those computed in Re . 

vibrational temperature. The properties oi unc f . . , b v an 

at 100 km. In comparison with the BSUV-2 data, the OH number density shown . S ' 

. . d the translational temperature shown in Fig. 6b is lower by a factor of a ou 

“of these 71 wil.'signiflcantly reduce the OH emission generated by ^ 77^7 

rr: r oh ^ was - - o. 3^, - ^—* -— 

predicted here for the cold-flow case would no. be visible from the Mau. facta,. 

For completeness, the internal mode temperatures of NH are shown in Fig 6c. The number ^ 
NH is even lower than that for OH while the translational temperatures are simi ar. 

that the NH emission would not be detectable from the ground. ^ 

The calculations did no. resolve the vibrationally hot NH, molecules that are formed y 
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the flow with ambient He or O. Also, it is estimated that there axe no vibrationally excited NH 3 particles 
leaving the nozzle. Nevertheless we can estimate the magnitude of the radiation from the cold gas arc-jet 
flow. The radiance estimates will be compared with those observed by NASA Infrared Telescope Facility 17 
of ammonia emission lines from the stratosphere of Jupiter. 

Based on the DSMC calculations the main source of radiation for the cold gas mode is if the ammonia 
itself radiates. The two IR radiation transitions that are considered are those at 2.9 and 10.5 fj, with emission 
rates of 2.2 and 15 transitions/sec, respectively. The radiance from a column of ammonia molecules may be 
approximated as, 


I = hc/\exp(-E/kT)ANAl 

where A is expressed in transitions/sec, N is the number density of the ground state of ammonia and from 
Fig. 5a is approximately 1 x 10 15 cm -3 , and A l is the layer thickness and is estimated from Fig. 5a to be 
20 cm. Figure 5b shows that the translational temperature of ammonia is very cold and is really less than 
25 K. However, assuming as an upper limit to radiation a higher value of 75 K, then the radiance of the gas 
at 10.5 (j. is about 3 x 10~ n W/cm 2 sr. This is a very small number and most likely is not detectable by a 
remote measurement. 

An example of a remote measurement of a similar infrared transition studied radiation from stratospheric 
ammonia on Jupiter. The data analyses of Kostiuk etal. 17 give a column abundance of ammonia of 2 x 
10 16 /cm 2 which is about the same as in the ammonia arcjet plume. However, the main difference is that 
the temperature of the Jovian ammonia is 150 K. Using these values in the expression above, a radiance of 
about 2 x 10~ 7 W/cm 2 sr is obtained. This value is four orders of magnitude larger than that estimated for 
the cold arcjet plume, thereby making it unlikely that a remote measurement could be made. The blackbody 
radiance at 150 K for this transition (which qualitatively represents the maximum radiance) in the same 
units is also of order 10 -7 W/cm 2 sr. Hence the estimate of the Jovian radiance is qualitatively correct. 

The plume flow calculations provide an estimate of the number density of vibrationallyhot NH 3 to be 
about-10® /cm 3 with temperatures in excess of 1,000 K. However, use of this estimate does not change the 
above conclusion that the plume radiance is too low to observe from a remote measurement. 

Conclusions 

The direct simulation Monte Carlo method was employed to model the nozzle, forward plume, and 
back flow plume regions of a cold flow ammonia arcjet. The numerical studies were performed for the flight 
conditions of the ESEX space experiment. 

The flow rate and nozzle geometry indicated near continuum conditions for the nozzle flow and this was 
confirmed by the DSMC computation. A significant boundary layer was formed in the nozzle, but the flow 
was in thermal equilibrium at the nozzle exit. 

The backflow region was well-resolved using a very large DSMC computation. Over 11 million particles 
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were employed In a simulation that ran in par* on 32 processors of an IBM SP-2. Significant fiuxe. of 

ammonia were predicted along all surfaces of the ESEX payload- He of 833 km was also sim- 

Interaction of the ammonia plume ^ ^'“^“^rgy to dissociate ammonia with the main 

ulated. The ambient ^ ^ ^ ultra . viol e, emissions from these molecules, 

products being N 2 , » the conclusion that none of the signals could be 

and the infra-red emission from the cold amm , ^ ^ ^ emi ssion ^ ^ ^ ^ 

detected-by an Earth-based, remote measurement. Shock Ultra Violet-2 reentry 

significantly smaller in magnitude compared to that detected in situ by the Bow Shock Ultra Violet 

flight experiment. 
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Rate coefficient, 
m 3 molecule -1 s -1 

k=AT'V -B /* r > 


Reaction ____ 

NH 3 + {He,0} -> NH 2 + H + {He,0} 
NH 3 + {He,0} -» NH + H 2 + {He,0} 
NH 3 + O -» NH 2 + OH 
NH 2 + {He,0} -» NH + H + {He,0} 
NH + { He,0 } -* N + H + {He,0} 

H 2 4- {He,0} H + H + {He,0} 

OH + {He,0} -*• O + H + {He,0} 

NH 3 + H NH 2 + H 2 

NH 2 + H <-» NH + H 2 

NH + H «-» H 2 + N 

NH 2 + NH 2 «-» NH 3 + NH 


1.41x10" 14 

9.13xl0 -15 

6.6xl0 -17 

9.13xl0 -16 

5.30xl0 -18 

5.50xl0 -18 

1.30xl0 -15 

2.10xl0 -16 

3.00xl0 -17 

1.66xl0 -18 

6.64xl0 -18 


E/fc, K 

46.507.3 
45,754.0 
4,455.0 

37.556.6 

34.826.6 

52.941.3 
50,978.0 
11,065.0 
2,626.8 
995.0 
2,800.0 


Reference 


Table 2. Mass flux in kg/m 2 /s at the measurement locations. 


Arcjet Mode 
Cold (present) 
Ignited (Ref. 1) 


QCM 1 

5.6xl0 -7 

3.0xl0 -8 


QCM 2 QCM 3 
3.5xl0 -1 ° 1.8xl0 -9 
1.2xl0 -9 1.2xl0 -9 


QCM 4 

2.5xl0 -8 

1.6xl0 -1 


Witness plate 
(with 8 solar cells) 

Radiometer (4 ptaces)^^^ 

Electromagnetic interference 
(EMI) antenna ^ 

Boom-mounted EMI antenna * 


EMI electronics 
Microbalance electronics 


Diagnostic package/ 
experiment platform 

Video camera and 
support structure 

Ouartz crystal 
' microbalances (4 places) 


tank 


Structure assembly 


Power 

conditioning unit 



Power integration unit^ 
Accelerometer and amplifier^ 
Camera electronics unit 


Propellant feed 
subsystem platform 


Command control unit 


Fig. 1. Schematic diagram of the ESEX'experiment. 
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Fig. 2a. Number density of ammonia in the nozzle exit. 



Fig. 2b. Temperatures of ammonia in the nozzle exit. 
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Fig. 2c. Velocity of ammonia in the nozzle exit. 
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Fig. 3a. Grid employed for the contamination simulation. 
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Fig. 3b. Number density contours for the contamination simulation. 
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Fig. 3c. Translational temperature contours for the contamination simulation. 
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Fig. 4a. Mass flu* profile along the surface of the ESEX payload starting at the nestle lip. 



Fig. 4b. Pressure profile along the surface of the ESEX payload starting at the noszle lip. 
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Fig. 6a. Number density of selected species along the plume centerline. 
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